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ABSTRACT 

The bispectrum provides a characteristic signature of gravitational instability that can be used 
to probe the Gaussianity of the initial conditions and the bias of the galaxy distribution. We 
study how this signature is affected by redshift distortions using perturbation theory and high- 
resolution numerical simulations. We obtain perturbative results for the multipole expansion of 
the redshift-space bispectrum which provide a natural way to break the degeneracy between bias 
and Q present in measurements of the redshift-space power spectrum. We propose a phenomeno¬ 
logical model that incorporates the perturbative results and also describes the bispectrum in the 
transition to the non-linear regime. We stress the importance of non-linear effects and show that 
inaccurate treatment of these can lead to significant discrepancies in the determination of bias 
from galaxy redshift surveys. At small scales we find that the bispectrum monopole exhibits a 
strong configuration dependence that reflects the velocity dispersion of clusters. Therefore, the 
hierarchical model for the three-point function does not hold in redshift-space. 


Subject headings: large-scale structure of universe; methods: numerical; methods: statistical 
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1. Introduction 

Understanding the origin of the large-scale distribution of galaxies in the Universe remains a fundamental 
goal of cosmology. Galaxy redshift surveys provide a unique resource for this endeavor. Recession velocities, 
however, are not exact indicators of galaxy radial distances, since gravitational instability induces peculiar 
velocities which contribute to galaxy redshifts. This departure from the Hubble flow causes the observed 
clustering distribution in redshift surveys to be statistically anisotropic. On large scales, the clustering 
pattern is compressed along the line-of-sight (the “squashing effect”) by coherent infall onto galaxy clusters 
and superclusters and outflow from voids and other underdense regions (Sargent & Turner 1977). On 
small scales, the velocity dispersion of virialized clusters creates the so-called “fingers of god” (FOG) by 
stretching out these structures along the line-of-sight (Jackson 1972). The anisotropy of galaxy redshift- 
space clustering provides valuable information about the galaxy velocity field and hence the underlying mass 
density field. In particular, measurements of redshift-space distortions can be used to determine the matter 
density parameter 

Modelling the redshift distortions of clustering statistics is a prerequisite for extracting information from 
redshift surveys (see Hamilton 1997 for a review). The effects of redshift distortions are well understood in 
linear perturbation theory: galaxy power spectrum measurements on large scales can be used to determine 
the combination (3 = where b is the linear bias factor (Kaiser 1987, Hamilton 1992). This degeneracy 

between bias and Q can be lifted in principle by using higher-order statistics, such as the bispectrum, to 
independently constrain b (Fry & Gaztahaga 1993; Fry 1994; Gaztahaga & Frieman 1994; Matarrese, Verde 
& Heavens 1997). However, to apply this technique to redshift surveys, redshift distortions of the bispectrum 
must be taken into account as well. Hivon et al. (1995) used second-order Lagrangian perturbation theory to 
calculate redshift distortions of the bispectrum monopole, and found them to be rather small, at most 10%. 
In this paper, we compute the redshift distortions in Eulerian perturbation theory and extend their results 
to the bispectrum quadrupole. The latter contains valuable information about fl, which complements that 
encoded in the power spectrum distortions. 

As has become evident in recent years from comparisons with numerical simulations, leading-order 
perturbative calculations—such as linear theory for the power spectrum—break down at larger scales in 
redshift space than in real space. In this work, we find a similar result for the bispectrum. Moreover, mildly 
non-linear effects are inevitably present in observations and must be adequately taken into account to take 
full advantage of the data (Cole, Fisher & Weinberg 1994, 1995; Hatton & Cole 1998). In particular, the 
next generation of galaxy redshift surveys, such as the Two Degree Field and Sloan Digital Sky Surveys, 
will provide their most precise measurements of redshift distortions on intermediate scales, where non-linear 
effects cannot be neglected. In order to improve upon the perturbative description of the redshift-space 
bispectrum, we develop a simple phenomenological model that extends the perturbative results to smaller 
scales, analogous to what has been done for the power spectrum multipoles. We compare our results with 
numerical simulations and find very good agreement even well into the non-linear regime. We find that the 
reduced redshift-space bispectrum Qs maintains its configuration dependence at all scales (in contrast to 
what happens in real space), due to anisotropies arising from the velocity dispersion along the line-of-sight. 
This effect, if not included, would otherwise compromise the program to use the configuration dependence 
of the bispectrum to constrain the bias from measurements in redshift surveys. 

The rest of this paper is organized as follows. In Section 2, we derive a fully non-linear expression 
for the density field in redshift space, which is then used to obtain the perturbative expansion order by 
order. Section 3 applies these results to yield the monopole and quadrupole moments of the redshift-space 
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bispectrum for different configurations, including non-linear biasing. The extension of these results into the 
non-linear regime via a phenomenological model is the subject of Section 4, where comparison with numerical 
simulations is made. Section 5 provides a hnal discussion and conclusions. 


2. Redshift-Space Mapping and Perturbation Theory 

In redshift space, the radial coordinate s of a galaxy is given by its observed radial velocity, a combination 
of its Hubble flow plus “distortions” due to peculiar velocities. The mapping from real-space position x to 
redshift space is given by: 

s = x-fu^{x)z, (1) 

where /(fl) ~ ® is the logarithmic growth rate of linear perturbations, and u{x) = —v{x)/{Ti.f), where 

v(x) is the peculiar velocity held, and 7Y(t) = {l/a){da/dT) = Ha is the conformal Hubble parameter (with 
FRW scale factor a(r) and conformal time r). In Eq. (I), we have assumed the “plane-parallel” (or inhnitely 
remote observer) approximation, so that the line-of-sight is taken as a hxed direction, denoted by z. 

The density held in redshift space, Ss{s), is obtained from the real-space density held S{x) by requiring 
that the redshift-space mapping conserves the number of galaxies, i.e. 


il + 6s)d^s = {l + 6)d^x . (2) 

Using the fact that d^s = J{x)d^x, where J{x) = 1 — fVzUz{x) is the exact Jacobian of the mapping in the 
plane-parallel approximation, we obtain: 


5{x) + l-J{x) 5{x) + zUz{x) 

=- W) -' 

The zeros of the Jacobian describe caustics in redshift space, the locus of points where the density held 
is apparently inhnite (McGill 1990). This surface is characterized in real space by those points which are 
undergoing turn-around in the gravitational collapse process, so their peculiar velocities exactly cancel the 
differential Hubble how. In practice, caustics are smoothed out by subclustering. 

We are interested in the Fourier transform of the density held contrast in redshift space. 


5s{k) = 


,k-s 

(27r)3 


Ss{s) 


d^x 


ik-x^ifk^Uz(x) 


5{x) + fVzUzix) 


( 4 ) 


This equation describes the fully non-linear density held in redshift space in the plane-parallel approximation 
and is the starting point for the perturbative approach. The term in square brackets describes the “squashing 
effect”, i.e., the increase in the clustering amplitude due to infall, whereas the exponential factor encodes the 
FOG effect, which erases power due to velocity dispersion along the line-of-sight. To obtain the perturbative 
expansion, we expand the second exponential in power series. 


Ssik) = ^ 

n—1 


d^kl . . . d^kn[6B]n S{ki) + ffil9{ki) 


(/M)"-3 

(n - 1)! k 2 


9{k2).. 



( 5 ) 


where [Join = ^(fc — fci — • • ■ — fen), with the Dirac delta distribution, the velocity divergence 

9{x) = V • u{x), and = ki ■ z/ki is the cosine of the angle between the line-of-sight and the wavevector. 
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In linear perturbation theory, only the n = 1 term survives, and we recover the well known formula due to 
Kaiser (1987) 

= (5(fe)(l +( 6 ) 

Equation (||) can be used to obtain the redshift-space density field beyond linear theory. For this purpose, 
we recall that the perturbative expansion for the un-redshifted density and velocity-divergence fields can be 
written (Fry 1984, Goroff et al. 1986) 


00 p 

6 (fe,r) = ^74-(r) / d^k,... 

n^l 


. .,fc„)(5i(fei)" 

■■Sl{kn), 

(7) 

CX5 p 

e{k,T) = Y,DUT) / d^k,... 

n=l 

J d^kn[6B]nGn{ki,. 

■■,kn)Si{ki)- 

■■Siikn), 

(8) 

where Di{t) is the density perturbation growth factor in linear theory. Here we have assumed that the 

nth- 


order growth factor Dn oc D^, which is an excellent approximation (see Scoccimarro et al. 1998, Appendix 
B.3). The perturbative kernels Fn and Gn are dimensionless, symmetric, scalar functions of the wave vectors 
{fei,..., fc„} which can be recursively found from the equations of motion. Similarly, in redshift space we 
can write the density field amplitude in the form 


OO 

6s{k,T) = Y,DUr) 

n—1 

To proceed, we assume a local, non-linear biasing scheme, in which the galaxy density contrast Sg is 
given as a Taylor series expansion in the underlying dark matter fluctuations <5 (Fry & Gaztahaga 1993), 


d kl . . . Id fcyi (fcl , . . . , fc^) (fcl) • • • (/Cyj). 


(9) 


Sg=J2^S-^^b{S + j6y2 + ...), ( 10 ) 

ml 

m—O 

where 7 = 62 /^, and b = bi is the linear bias. If the galaxy distribution is unbiased, then 6=1 and bm = 0 
for TO > 1. The term 60 ensures that (Sg) = 0; since it does not enter into connected correlations, we will 
neglect it in what follows, as assumed in the last expression in Eq. We note that the model in Eq. ®) 
is a rather simple bias prescription, and one could generalize the results below to more complex biasing 
schemes (e.g., non-local, time-dependent, stochastic, etc). 

From Eqs. (|) -(p^, we can now read off the redshift-space kernels for the galaxy density field. In 
particular the linear, second- and third-order perturbation theory (PT) kernels are 


= (& + V), 

= bF2{ki,k2) + fn'^G2{ki,k2) + ^-1^ + /mI) + ^(^ +/Mi) +17’ 

2 L K 2 -I 2 

= bF^^\ki,k2, fca) + fti^Gi^\ki,k2, fea) + ffik\bF^^\k,,k2) + ^ 2 ) 


( 11 ) 

( 12 ) 

in 

ks 


+ffik{b + Jy,l)I^Gi^\k2M) + + 362Fi^^(fci, 

/C 23 2 /C 2 AC 3 


fc.)+ 1 , 


Z,(k} 

Z2{ki,k2) 

Z3(fci,fe2,fe3) 


2 


( 13 ) 







- 5 - 


where we denote fi = k ■ z/k, with fc = fci + ... + A;„, and = ki ■ z/ki. As above, F 2 and G 2 (given 
below) denote the second-order kernels for the real-space density and velocity-divergence 


in Eqs. (22 


fields, and similarly for F 3 and G 3 . Note that the third order kernel must still be symmetrized over its 
arguments. One can similarly obtain the PT kernels in redshift space to arbitrary higher order. 


We note that there are two approximations involved in this procedure: one is the mathematical step of 
going from Eq. (Q) to Eq. (^, which approximates the redshift-space mapping with a power series; the other is 
the PT expansion itself (i.e., the expansion of 6{k) and 9{k) in terms of linear fluctuations Si{k)). Therefore, 
one is not guaranteed that the resulting PT in redshift space will work over the same range of scales as in 
real space. In fact, we will find that, in general, PT in redshift space breaks down at larger scales than in real 
space, because the redshift-space mapping is only treated approximately, and it breaks down at larger scales 
than does the perturbative dynamics. In particular, a calculation of the one-loop power spectrum in redshift 
space using Eqs. does not give satisfactory results. To extend the leading-order calculations such as 

those in this work, one must treat the redshift-space mapping exactly and only approximate the dynamics 
using PT (Scoccimarro, Couchman & Frieman 1998; hereafter SCF). 

The calculation of redshift-space statistics in Fourier space proceeds along the same lines as in the 
un-redshifted case. To leading (linear) order, the redshift-space galaxy power spectrum reads (Kaiser 1987) 


Ps{k) = Pg{k) (1 + /3m') 2 = ^ a, Pii^^) Pgik), (14) 

where Pg{k) = b^P{k) is the real-space galaxy power spectrum, P{k) is the linear mass power spectrum, and 
jd = f /b K, jb. Here Pi{pl) denotes the Legendre polynomial of order £, and the multipole coefficients 
are (Hamilton 1992; Cole, Fisher & Weinberg 1994) 

oq = 1 +-jd +, a2 = -P+-P'^, 04 = —(d^; (15) 

3 5 3 7 35 

all other multipoles vanish. Equation (0) is the standard tool for measuring ft from redshift distortions 
of the power spectrum in the linear regime; in particular, the quadrupole-to-monopole ratio Rp = 02/03 
should be a constant, independent of wavevector k, as k =>■ 0. Note, however, that in these expressions H 
appears only through the parameter fd, so there is a degeneracy between n and the linear bias factor b (it is 
impossible to distinguish an unbiased low-density Universe from a biased high-density model). 

Given the second-order PT kernel in redshift-space, the leading-order (tree-level) galaxy bispectrum in 
redshift-space reads 


Bs{ki,k 2 ,kz) = 2 Z 2 {ki,k 2 ) Zi{ki) Zi{k 2 ) P(fci) P{k 2 ) -bcyc., (16) 

which can be normalized by the power spectrum monopole to give the hierarchical three-point amplitude in 
redshift space, Qs, 


Qs{ki,k2,k'i) = 


Bsiki,k2,k3) 


Bs{ki,k2,k3) 


(17) 


Ps{ki) Ps{k2) + eye. al [Pg{ki) Pg{k2) -f eye.) ’ 

where “eye.” denotes a sum over permutations of {ki,k 2 ,k 3 \. Note that the hierarchical amplitude Qs 
is independent of power spectrum normalization to leading order in PT. Since, to leading order, Qs is a 
function of triangle configuration which separately depends on Cl, b, and 62 , it allows one in principle to 
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break the degeneracy between and b present in measurement of the power spectrum multipoles in redshift 
space. 


3. Results: The Redshift-Space Bispectrum in Perturbation Theory 

Due to redshift distortions, the bispectrum becomes a function of five variables: three of them describe 
the shape of the triangle (e.g., the sides fci, k 2 , and the angle 9 between them, with cos0 = fci • ^ 2 ), and the 
two remaining variables characterize the orientation of the triangle with respect to the line-of-sight, which 
we take to be the polar angle to = cos“^ /r of fci, and the azimuthal angle (j) about fei: 


/il = /i = cosw = fci 




= cos 9 — \/{l — /r^) sin 9 cos (j), 


k3 = ^9-2- (18) 


The dependence introduced by redshift distortions can be conveniently described by decomposing the 

tree-level bispectrum Bg in Eq. (0) into spherical harmonics, 


00 i 

Bs{k,,k2,k3)=J2 ( 19 ) 

l—O ra——t 

Bispectrum multipoles are non-vanishing for even £, up to £ = 8, m = 6. We now explore the 

information encoded in these multipoles for different triangle configurations. 


3.1. Non-Equilateral Configurations 

Rather than working with the full multipole decomposition (€, 771 ), here we concentrate for simplicity 
on the m = 0 multipoles. In this case, which corresponds to averaging over </), it is more convenient to 
decompose the bispectrum in Legendre polynomials 


R,(fci,fc2,fc3) = ^ Bi^\ki,k2,9) PtifJ.). (20) 

e=o 

We can write the different contributions as 

= Pg{ki)Pg{k2) ^E2(fci, fc2) 2 ?sQ 1 + G'2(fcl, fc2) ^SQ2 + ^NLB + ^FOg) + ^JC., (21) 


where the second-order real-space kernels are given by 


F2(fcl,fc2) 


G2(fcl,fc2) 


5 X f ki fc2 \ 

7^2\Y2^YJ 

3 X 

7^ 2 VA:2 kj 


2x^ 

~T' 

4x^ 


( 22 ) 

(23) 


with X = cos 9. The first term in Eq. ( pT[ ) represents the effect of the linear squashing (SQl), the next 
term describes second-order squashing (SQ2), the third contribution is due to non-linear biasing (NLB), 
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and the last term describes the effect of damping due to velocity dispersion, the FOG effect. The different 
contributions to the monopole {£ = 0) and quadrupole {£ = 2) are given by (recall that /? = //fo, 7 = 62 /^) 

T>e=o — ^ + 10 /3 + + 2 x'^) 

= 1^6 
£^2 _ 2/3 (7 + /3 + 21 + 11 /3a:^) 

^SQl = 

= 2f3{35ki^ +28pki^ + 3(3“^ ki^ + 35 + 28/3 + 3 /S^ + 70 kik2X + 84(3kik2X + 

18/3^ fci a; + 14/3 ki^ x^ + 12 (3^ ki^ x'^ + I4j3k2^ x^ + 12 /3^ k 2 ^ x^ + 12/3^ fci ^2 
/(105fc3^6) (26) 


(24) 

(25) 


= 2/3(14/72 + 13/3^1^ +/32A:i 2 -7fc2^ -5/3fc2^ + 28A:ifc2a; + 30/3A:iA:2a; + 6/32fciA:2a; + 

11 /3/ci 2 a;2 + 9/32 fci^ a;2 + 21 ^ 2 ^ a;2 + 29/3/c2^ a;2 + 6/32 ^ 2 ^ a;2 + 18/3 fci ^2 x^ + 14/32 fci ^2 x^ + 
4(3^ k2^x^)/(21 ks^b) (27) 


P 


t=o 

NLB 


V 


£=2 

NLB 


7 (15 + 10/3 +/32 + 2/32 x2) 
15b 

(3j (7 + /3 + 21x 2+ 11/3x2) 
216 


(28) 

(29) 


= /3 (210 fci fc2 + 210 /3 fci fc2 + 54 /32 fci fc2 + 6 ki *2 + 105 /c? x + 189 /3 /c? x + 99 (3^klx + 
15/3^ fc2x + 105 fc2x + 189/3/c2a; +99/32 fc2^_pi5^3 + 168/3/7 *2 a;2 + 

216 /32 /ci /c2 x2 + 48 /3^ /ci /c2 a;2 + 36 0^ kf x^ + 20 /3^ kl x^ + 36 (3'^ /c| x^ + 20 /3^ /c| x^ + 

16/3^/ci/c2 a;‘*)/(315/ci/c2) (30) 


V00Q = (i (231 /7 k 2 + 330 /3 /ci ^2 + 99 /32 /ci ^2 + 12 /33 /ci /c 2 + 462 /72 x + 891 /3 372 x + 528 0^ k0 x + 

75 (3^ k0 X + 462 x + 594 /3 /c 2 ^ a; + +198 0^ /c 2 ^ x + 30 /33 /c 2 ^ a; + 693 /ci k 2 X^ + + 

2046 (3 ki k2 x"^ + 1485 0^ ki /c2 a;2276 (3^ ki /c2 ^2 + 297 (3 k0 x^ + 462 (3^ k0 x^ + 205 (3^ k0 x^ + 
594 /3 /C2^ a:^ + 792 0^ k2^ x^ + 190 /3^ /c2^ a;^ + 396 ki /c2 a:'‘ + 272 /3^ /ci /c2 a:'‘ + 

60P^ k2^ x01(693 kik2) (31) 

The result for the monopole without bias, when 6=1 and 7 = 0 , agrees with that given in Hivon et al 
(1995), where they also include the (extremely weak) fl-dependence due to second-order PT growth factors 
(i.e., the small deviation from D 2 oc Di). This agreement may seem surprising, since they claim to include the 
radial character of redshift distortions, whereas our calculations are done in the plane-parallel approximation. 
In fact, although intermediate expressions in Hivon et al (1995) do preserve the radial character of redshift 



distortions in the large-volume approximation [equivalent to replacing z by r in Eq. (Q)], the angle averaging 
is done instead assuming that r is a fixed direction, that is, in the plane-parallel approximation. Our results 
generalize theirs to include biasing; more importantly, the result for the bispectrum quadrupole is new. 

Figure 1 illustrates these results. The top left panel shows the amplification of the bispectrum due to 
redshift-space distortions, showing the ratio /B, for ^ = 0 and £ = 2, as a function of angle 9, 

for configurations with ki/k2 = 2 and scale-free initial conditions, P{k) (x k'^ {n = —2 solid curves, and 
n = 0 dotted curves). The monopole shows a small (< 20%) change of the bispectrum due to redshift 
distortions, as first shown by Hivon et al. (1995). The quadrupole, however, shows a higher amplification 
for co-linear configurations (near 0 = 0 , 7 r, where fci||fc 2 ) and a suppression at equilateral-like configurations. 
This is to be expected, since the quadrupole is generated solely by redshift distortions, and the dominant 
peculiar velocities in PT are co-linear with density gradients. The top right panel shows the bispectrum 
quadrupole-to-monopole ratio, Rb = Bf'^/B^^^ for r = ki/k 2 = 10 , 2,1 (top to bottom), in fl = 1 (top set) 
and n = 0.3 (bottom set) models, for n = —2 (solid) and n = 0 (dotted). As in the power spectrum case, 
Rb increases with 17, and it shows the expected configuration dependence; it is relatively insensitive to the 
shape of the power spectrum. 


The bottom right plot in Figure 1 quantifies how redshift distortions alter the relation between galaxy 
and matter correlations. For the local, non-linear bias model of Eq. (0), the 3-point hierarchical amplitude 
for galaxies satisfies Qg = {Q + j)/b in the absence of redshift distortions (Fry & Gaztahaga 1993, Fry 1994). 
As can be seen from the different terms in Eq. (21), this simple relation no longer holds in redshift space. 
The solid curves show the exact PT result in redshift space, while the dotted curves correspond to (wrongly) 
assuming that Qg s = {Qs + 'y)/b holds in redshift space, i.e., that the operations of bias and redshift-space 
mapping commute. Here we have set the linear bias to 5 = 2, the non-linear bias 7 = 1/2,0, —1/2 (from top 
to bottom), and show configurations with ki/k 2 = 2, for a model with linear spectral index n = —2. We see 
that the assumption that bias and redshift distortions commute is accurate at about the 20 % level. 


3.2. Equilateral Configurations 

For equilateral triangle configurations (fci =^2 = ^ 3 ), the leading-order PT expression for the redshift- 
space bispectrum, Bg eq, is 


Bg eq(/i, ly ) = [ Pg { k )]‘^ (256 -f 448 7 -f 288 ( 3 + 22 Ab ( 3 + 448 /37 + 72 /32 ^4 168 6/ + 

84/327/ -h4/3V® - 7 6/3"* / -f 288/31/2 -H 2245/3 j/2 -H 448/371/2 _ 288/3/^2 _ 

224 5 /3 /r2 j/2 _|_ 02 02 ^2 _|_ 2^0 ^ ^2 02 ^2 _ ^ ^ ^2 ^2 _|_ 2^0g 02 ^ 02 ^2 _ ^2 ^4 ^2 _ 

33 6 5 /32 / 1/2 — 24 /3^ jy2 _ 2^0g ^2 ^ / j^2 _|_ 24 ^3 ^6 ^^2 gg ^ ^4 jy2 _ gg g / ^2 _|_ 

72 /321/4 -H 168 5 /321/4 + 84 144 ^2 _ ggg g ^2 ^2 ^ ^ gg ^3 ^2 ^ _ 

168 332 7 7(2 1/4 + 72 p'^ 7(4 0 _|_ 408 5 02 ^4 0 — 'j2 p^ 0^ 0 — 21b /3'^ 0^ j/4 -1- 84 /32 7 0^ 0 + 

36 P^ 74® 1/4 -I- 42 5 /3^ 7(® i/"^ — 215 P'^ 0 — 63 b /3^ 0^ 0 + 189 5 p '^ 0^ 0 — 

189 5/3^7(®(/®-f 63 5/3V®i^®) x 3/(448 5) (32) 

where z/ = cos^. For comparison, in real space, the PT bispectrum is simply Heq = (12/7)[Pg(fc)]2 (Fry 
1984). The bottom left panel in Fig. shows the ratio of the redshift- and real-space equilateral bispectra in 
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leading order PT, Ab = Bg For /r = 0 and (/> = 7 r/ 2 , the triangle lies in the plane perpendicular 

to the line-of-sight, and there are no redshift distortions {Ab = !)• The maximum distortion occurs for /r = 0, 
(p = TT, with an enhancement factor Ab = 189/48 = 3.94. A similar enhancement {Ab = 1005/256 = 3.93) 
is obtained for the configuration with fi = 1, independent of tp. 

Averaging over p, the equilateral redshift-space bispectrum becomes 


Bseq{f^) = [Pg(fc)]2 (^4096 + 2304/3+ 1792 6/3 + 432/32 + 1008 6/3^ + 7168 7 + 3584/37 + 504/32 7 + 

2304 /3 /r2 + 1792 6 /3 /r2 + 288 + 67 2 6 /32 ii2 + 216 /3V^ - 315 6 /3^ + 3584 Pjf / + 

336 /32 7 ^2 432 ^2 ^4 4QQg _ 534 ^4 349 5 ^4 ^4 5Q4 0 ^ ^ 

472/3V®-413 6/3V® - 203 6/3^/) x 3/(7168 6) (33) 

Decomposing into Legendre polynomials, Bg eq {^^) = EZo Pi{fJ.), and defining the equilateral hierar¬ 

chical amplitude in redshift-space, Qi^eZ = Ss^eq*V3[Ts(fc)]^, we obtain 



5 (2520 + 4410 7 + 1890 /3 + 2940 7 /3 + 378 /32 + 441 + 9/3^ + 1470 6 /3 + 882 6 /32 

98 6(15 + 10/3 + 3/32)2 


14 6/3^) 


(34) 

This result shows once again that in redshift space, Qg^g 0 {Qg + 7)/6 (although, as shown in Fig. 1, it is 
not a bad approximation). In the absence of bias (6 = 1, 7 = 0), Eq. (^) yields 


Q 


{£= 0 ) 

s eq 


5 (2520 + 3360/ + 1260 p + 9f-Uf) 
98(15 + 10/+ 3/2)^ 


(35) 


which approaches the real-space result Qeq = 4/7 = 0.57 in the limit / ~ ® ^ 0. On the other hand, for 

/ = D = 1 , we have Qi°eq = 0.464: for these configurations, the hierarchical 3-point amplitude is suppressed 
by redshift distortions. For the quadrupole-to-monopole ratio of Bgq we have 


eq — 


5 (4158 /3 + 3234 7 /3 + 1188 /32 + 693 j/3^ + 33/3^ + 32 3 4 6 /3 + 27 7 2 6 /32 - 56 6 (3*) 

22(2520 + 4410 7 + 1890 /3 + 2940 7 /3 + 378 /32 + 441 7 /32 + 9 /33 + 147 0 6 /3 + 88 2 6 /32 - 14 6 (3^)' 

(36) 

provide two constraints upon 


For no bias and 0 = 1, this gives Rb eq = 11329/31394 = 0.36. Eqs. (34 


{b,(3,0 that are independent of the power spectrum. Thus, these relations for eZ and i?B eq, together 
with the power spectrum quadrupole-to-monopole ratio Rp from Eq. (H), could in principle be inverted to 
obtain O, 6 , and 62 at large scales, independent of the initial conditions. In practice, Rp will determine /3, 
while Q^seZ and Rb eq then constrain the 6,7 parameter space. Moreover, Rp is very insensitive to 6 and 
only mildly sensitive to 7 , so that there will typically be some residual degeneracy between 6 and 7 . To break 
this degeneracy, one could use the shape-dependence of B, i.e., non-equilateral configurations (Fry 1994). 


4. Results: The Redshift-Space Bispectrum in Numerical Simulations 
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4.1. Numerical Simulations 

The numerical simulations used in this work correspond to cluster-normalized cold dark matter (CDM) 
models run by the Virgo collaboration (see e.g., Jenkins et al. 1998 for more details). We present results 
for standard CDM (hereafter SCDM, with = 1, h = 0.5) and “Lambda” CDM (hereafter ACDM, with 
n = 0.3, Ua = 0.7, h = 0.7), although we have also done measurements in rCDM (0 = 1 with power 
spectrum shape parameter V = 0.21) and open CDM (O = 0.3, h = 0.7) models. We found the results 
for SCDM and ACDM to be representative of the full set of simulations. These simulations each contain 
256^ particles in a box 240h~^ Mpc on a side. They were run with an adaptive P^M code (Couchman, 
Thomas & Pearce, 1995; Pearce & Couchman, 1997), and the initial conditions were set using the Zel’dovich 
approximation on a “glass” (see, e.g. White 1996). We have also made measurements in a set of smaller 
simulations run by the Hydra Consortium (Couchman, Thomas & Pearce, 1995), with very similar results. 

We use the plane-parallel approximation to construct the redshift-space maps, so that the periodic 
boundary conditions in the simulation box are preserved by the redshift-space mapping. We can therefore 
continue to use the FFT to obtain the density field in Fourier space. This also has the advantage that our 
N-body measurements are done using the same approximation as our perturbative calculations. We consider 
four different observer’s directions, in all cases misaligned with respect to the simulation box sides to avoid 
spurious effects. Error bars in the plots denote the dispersion in the values obtained for these four observers. 
To measure the bispectrum in the numerical simulations, we use the method described in Scoccimarro et 
al. (1998), developed by S. Colombi, with minor improvements such as using TSC instead CIC to do the 
interpolations. 


4.2. Phenomenological Description of Non-Linear Redshift Distortions 

In order to describe the non-linear behavior of the redshift-space bispectrum, we introduce a phenomeno¬ 
logical model to take into account the velocity dispersion effects and to quantify their importance. We follow 
the standard approach in the literature (Peacock & Dodds 1994), in which the non-linear distortions of the 
power spectrum in redshift-space are written in terms of the linear squashing factor and a suitable damping 
factor due to the pairwise-velocity distribution function 


Ps{k) 


Pgik) 


(1 + /3m^)^ 

[1 -I- (A:/x(t„)2/2]2 


(37) 


Here ay is a free parameter that characterizes the velocity dispersion along the line-of-sight. This Lorentzian 
form of the damping factor is motivated by empirical results showing an exponential one-particle velocity 
distribution function (Park et al. 1994); comparison with N-body simulations have shown it to be the best 
phenomenological model so far (Cole, Fisher & Weinberg 1995; Hatton & Cole 1998). It is straightforward to 
obtain the multipole moments of Ps{k) in this simple model (Cole, Fisher & Weinberg 1995). In comparing 
with numerical simulations, we use the quadrupole-to-monopole ratio statistic = 02/00 to fit dt,. 


We propose a similar ansatz to model non-linear distortions of the bispectrum: 


Hs(fci,fe2,fe3) = 




(38) 


[1 + [(fciMi)" + (fc2M2)2 + (A:3/X3)2]V2/2]" 
where iJf'^(A;i, fc 2 , ^ 3 ) is the tree-level redshift-space bispectrum derived in the last section. Note that we 
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have written the triplet velocity dispersion along the line-of-sight in terms of the pairwise velocity dispersion 
parameter tT„ and a constant a which reflects the configuration dependence of the triplet velocity dispersion. 
As noted above, is determined from simulations solely using the power spectrum ratio i?p; the parameter 
a is then fitted by comparison with the monopole-to-quadrupole ratio of the equilateral bispectrum, Rb, 
measured in the simulations. The simpler phenomenological model in which the supression factor in the 
denominator in Eq. (^) is replaced by 1 + + (^ 21 ^ 2 )^ + (fc3M3)^]o’^/2 does well in comparison with 

the N-body results for configurations where k 2 /ki = 2, but does not reproduce the results of the equilateral 
quadrupole-to-monopole ratio Rb- In what follows we find good fits to the measured bispectrum using 
a ~ 2 for equilateral configurations and a ~ 3 for fci/fc 2 = 2 configurations, independent of cosmology (at 
redshift z = 0). For z = 1, we find that a = 1 (SCDM) and a = 1.75 (ACDM) reproduce the shape of 
Rb- This is similar to the modelling of the equilateral three-point function by Matsubara (1994), who found 
it necessary to have a triplet velocity dispersion parameter ra l.Tay to fit the N-body simulations. We 
note here that we have not attempted to do any best-fit procedure to obtain (T„ and a: these values were 
simply obtained by inspecting by eye the results for Rp and the bispectrum. We also note that eqn. (pq) 


reduces to the perturbative expression at large scales (k 
[e.g. Eqs. 


0), so it has the correct limiting behavior 




4.3. Comparison with Numerical Simulations 


Figure g shows the results from the phenomenological redshift-distortion model and how it compares 
with N-body simulations. Error bars have been suppressed for clarity but can be guessed from the dispersion 
of the data points. For the SCDM model, from Rp we obtain = 6 (in units of Ho=100 h km/s/Mpc) for 
z = 0 and ay = 2 for z = 1. Using this result, we can predict the equilateral bispectrum ratio Rp in the 
phenomenological model by taking multipoles of Eq. (38). The results are shown as the solid curves in Fig. 2, 
showing excellent agreement with the N-body simulations. Similarly, for the ACDM model, from fitting Rp 
we obtain ay = 5.5 for z = 0 and cr^ = 4 for z = 1. The results for i?B are again in very good agreement 
with the numerical simulations. Note the similar behavior of Rp and Rp as a function of scale: the FOG 
effect, which creates a deficit of pairs along the line-of-sight, also causes a decrease of triplets compared to 
those in the perpendicular direction. In this sense, it is a reassuring check to note that in every case the 
zero-crossings of Rp and i?B are very close to each other in scale. 


The top left panel of Fig. 3 shows the power spectrum amplitude A(k) = 47 rk^F(k) in real (squares) 
and redshift (triangles) space for SCDM, with the corresponding linear theory predictions. In the figures, 
perturbative predictions are denoted by “PT” (real space) and “PTs” (redshift space). At large scales, the 
squashing effect boosts the power (Kaiser 1987), whereas at small scales the FOG effect erases power in the 
monopole of the power spectrum. The top right panel shows the equilateral hierarchical amplitude Q in real 
(squares) and redshift (triangles) space, with the corresponding predictions of tree-level PT given by Eq. (^) 
in redshift space and Qeq = ^jl in real space. Note the similarity to the power spectrum results: at small 
scales the FOG effect maintains the redshift-space equilateral amplitude close to its perturbative value, while 
in real space the amplitude grows by a large factor. This result, or the analog of this result for the skewness 
S 3 {R) (i.e., the fact that S 3 in redshift space is less scale-dependent than in real space), has led previous 
authors to conclude that higher-order correlations are ‘more hierarchical’ in redshift-space than in real space 
(Lahav et al. 1993; Matsubara 1994, Matsubara and Suto 1994; Suto and Matsubara 1994; Bonometto et al. 
1995; Ghigna et al. 1996). However, this interpretation is in fact not correct: the FOG effect at small scales 
decreases the equilateral amplitude with respect to the real space value and shifts power to the co-linear 





- 12 - 


configurations. Therefore a strong configuration dependence of the bispectrum is expected in redshift space, 
even when the real-space correlations are well described by the hierarchical model (Scoccimarro et al. 1998). 
This is exactly what we find in numerical simulations (bottom right panels in Figs. 3 and 4), and what the 
model in Eq. (|^ ) predicts. 

The bottom panels in Fig. 3 show a comparison of simulations to the predictions of PT and the phe¬ 
nomenological model of Eq. (|^ ) for configurations with k 2 /ki = 2, for two different scales: fci = 10 kf (left) 
and ki = 20 kf (right), where the fundamental wavenumber of the simulation box is kf = 27r/240 hMpc“^. 
In the bottom left panel, where the amplitude (A) ft! A(fc 2 ) = 2.5 (thus somewhat beyond the weakly non¬ 
linear regime), the tree-level PT calculation in real space (dashed curve) predicts the form of the N-body 
results for Q reasonably well, although not to high accuracy. However, the redshift-space PT counterpart 
(dotted curve) does not provide a good description of the simulation results for Qg- This illustrates the 
point that PT tends to break down on larger scales in redshift space than in real space. On the other hand, 
the phenomenological model of Eq. (|^ ) (solid curve) describes the N-body results for Qg very well. If one 
were instead to use the redshift-space PT prediction [Eq. (^] at these scales to probe the bias, one would 
mistakenly conclude that the bias is b = 1.7 instead of 1. 

A similar set of plots is shown in Figure 4 for the ACDM model, corresponding to ki/kf = 5,10,20. 
Even at the largest scale probed, corresponding to a wavelength A ft 50 Mpc//i, the tree-level PT prediction 
in redshift space is a poor fit to the simulation results; using it, one would erroneously predict an effective 
bias b = 1.4, a quite significant discrepancy. The bottom right panel nicely illustrates our point about 
correlations in the strongly non-linear regime: whereas Q is very close to hierarchical (nearly configuration- 
and scale-independent), the redshift-space amplitude Qg displays a strong configuration dependence. The 
phenomenological model (solid curve) does an excellent job in predicting Qg^ even at this stage of considerable 
non-linearity. 

In a recent paper, Verde et al. (1998) independently proposed a somewhat different phenomenological 
model to take into account redshift-space distortions of the bispectrum, Sg{k) = + {kfiav)"^ 

which makes their damping factor in Eq. ( |3^ the square root of ours. We find that the resulting power 
spectrum quadrupole-to-monopole ratio i?p does not have the correct shape when compared to numerical 
simulations, although it can be made to fit the low-fc behavior. If we fit ay to the simulation Rp using their 
model, we find that the predicted Qs does not provide a very accurate fit to the simulation results. In 
particular, comparing their model to the simulation results in Figs. 3 (bottom left) and 4 (top right) would 
lead one to erroneously conclude that there is an effective bias of 6 = 1.5 and b = 1.4, respectively. 

The Verde et al. phenomenological model can however be improved as follows (Heavens 1998, private 
communication). The power spectrum quadrupole-to-monopole ratio can be fixed by introducing a scale 
dependence of the velocity dispersion parameter ay —> av{k). Furthermore, using the non-linear power 
spectrum in the perturbative formula for the bispectrum, Eq. (p^, and removing the wavevectors close to 
the line of sight leads to Qg values that agree with numerical simulations. However, this statistic is not 
exactly equal to the monopole average since some configurations are not included in the averaging 

procedure. We conclude, therefore, that the phenomenological model of Eq. (|8|) is the simplest way to 
accurately extract bias parameters from redshift-space data. 
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5. Conclusions and Discussion 

We have derived the fully non-linear expression for the density field in redshift space in the plane-parallel 
approximation and applied this result to obtain the redshift-space bispectrum to leading order in PT. This 
derivation in Eulerian space is both simple and transparent and can in fact be extended beyond leading 
order to understand the large-scale FOG effect from a dynamical point of view (SCF). We computed the 
multipole expansion of the bispectrum in redshift space, in particular the monopole (Hivon et al. 1995) 
and quadrupole, including the effects of non-linear biasing. These perturbative results provide in principle a 
means of breaking the degeneracy between bias and O present in measurements of the redshift-space power 
spectrum on large scales. 

We found that, even at relatively large scales (e.g., 50 h“^Mpc) the leading-order PT predictions in 
redshift space are not accurate enough for detailed modelling and can lead to systematic errors in the deter¬ 
mination of the bias. This is similar to what happens with the power spectrum: redshift-space statistics are 
more strongly affected by non-linearities than their real-space counterparts. In this respect, it is interesting 
to note that even in the linear dynamics, the exponential factor in Eq. i) can lead to a FOG effect at 
large scales (SGF). Therefore, the long range of the FOG effect seen in numerical simulations, should not 
be attributed exclusively to virialized clusters. This behavior agrees with the results of Taylor & Hamilton 
(1996) and Fisher & Nusser (1996), who found that they could qualitatively reproduce the shape of Rp seen 
in numerical simulations by using the Zel’dovich approximation, although quantitatively its use can lead to 
systematic errors (Hatton & Gole 1998). The large-scale FOG effect has been recently connected to redshift- 
space shell-crossing due to coherent infall by Hui, Kofman & Shandarin (1998). All of these results strongly 
suggest the possibility of extending the leading-order PT results for the power spectrum and bispectrum to 
smaller scales by treating the redshift-space mapping in Eq. (jj) exactly and approximating the dynamics 
using PT (SGF). 

To correct for these non-linear effects, we have developed a phenomenological model that convolves the 
PT results with a damping factor due to velocity dispersion. Gomparison with high-resolution numerical 
simulations shows that this model, calibrated using the power spectrum distortions, works very well and 
provides a quantitative tool that can be applied to galaxy redshift surveys. 

In the strongly non-linear regime, the reduced bispectrum is very close to hierarchical (configuration- 
independent) in real space. However, in redshift space we found that there is a strong configuration de¬ 
pendence that reflects the velocity dispersion of clusters. This helps to understand the deviations from the 
hierarchical ansatz reported recently in the LGRS redshift-survey three-point function (Jing & Horner 1998). 
For equilateral configurations, on the other hand, the FOG effect suppresses the scale-dependence of Qs eq 
compared to that of Qeq, and this can be misinterpreted as a sign of hierarchical correlations in redshift 
space. A similar effect occurs with the skewness factor S 3 , which in redshift space shows a remarkable 
scale-independence from weakly non-linear scales into the non-linear regime (Hivon et al. 1995). Although 
Qeq and S 3 in redshift space remain very close to their PT values well into the non-linear regime, one must 
be careful using the standard non-linear biasing prescriptions beyond the regime of validity of tree-level PT, 
since non-commutativity of redshift distortions and biasing may be more complex than that exhibited at 
large scales. 
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Fig. 1.— Top left panel shows the redshift-space tree-level bispectrum multipoles, jB, as a 

function of angle 0 = cos“^ fci • k 2 , for configurations with ki/k 2 = 2 and scale-free initial spectra, P{k) oc 
fc" (n = —2 solid, n = 0 dotted). Top right panel shows the bispectrum quadrupole-to-monopole ratio 
Rb = for r = ki/k 2 = 10, 2,1 configurations as a function of 0 (n = —2 solid, n = 0 dotted) for 

n = 1 and n = 0.3. The bottom left plot shows the ratio of the redshift-space to real-space bispectrum 
for equilateral triangles, as a function oi fi = ki ■ z for different azimuthal angles (p. The bottom right plot 
shows the hierarchical amplitude Qs for ki/k 2 = 2 configurations with linear bias 6 = 2, and quadratic bias 
7 = b 2 /b = 1/2,0,—1/2 (top to bottom), for n = —2. Solid curves represent the PT result, and dotted 
curves assume that bias and redshift-space mapping commute. 
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Fig. 2.— The quadrupole-to-monopole ratio for the redshift-space power spectrum, i?p, and equilateral 
bispectrum, Rb- The symbols correspond to AP^M 256^ particle N-body simulations (see text) at redshift 
z = 0,1 averaged over four different observers, for SCDM (top) and ACDM (bottom). Squares denote 
measured values of i?p, and triangles correspond to values of the ratio i?B. The dashed and solid curves 
show the predictions of PT convolved with an exponential velocity dispersion model, Eqs. (|3^^^. 
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Fig. 3.— The top left panel shows the power spectrum amplitude, A(/c) = A'Kk^P{k) in real (squares) and 
redshift (triangles) space. The top right panel shows the equilateral hierarchical amplitude Qeq in real-space 
(squares) and the monopole of the redshift-space amplitude, Qs eq (triangles). The bottom panels show the 
hierarchical amplitude Q for k^jkx = 2 configurations, as a function of the angle 0 between fci and in real 
(squares) and redshift space (triangles), for two different scales. In the bottom left plot, weakly non-linear 
distortions decrease the configuration dependence of Qa^ similar to the effects of bias in leading order PT. 
At smaller scales (lower right panel), the configuration dependence of Qs is greatly enhanced, reflecting the 
velocity dispersion of virialized clusters. 
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Fig. 4.— Same as Fig ||, for the ACDM model. In the bottom right panel note the difference between 
the real- and redshift-space hierarchical amplitudes. Whereas Q in real-space is approximately constant, its 
redshift-space counterpart depends strongly on configuration. Therefore, the hierarchical ansatz does not 
hold in redshift space. 







































































